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O : ABSTRACT 

: 

We analyze low and high-latitude fast solar wind data from the Ulysses space- 

2 ■ craft from 1992 to 1994 using a a systematic method to analyse the anisotropic 

c/3 . content of the magnetic field fluctuations. We investigate all available frequen- 

ce ' a 

• • . cies, 1 — 10^ Hz, for both high and low-latitudes datasets and are able to quantify 



the relative importance of the anisotropic versus the isotropic fluctuations. We 
analyse, up to sixth order, longitudinal, transverse and mixed magnetic field cor- 
relations. Our results show that strongly intermittent and anisotropic events are 
present in the solar wind plasma at high frequencies/small scales, indicating the 
absence of a complete recovery of isotropy. Anisotropic scaling properties are 
compatible for high and low-latitude data, suggesting a universal behaviour in 
spite of the different rate of evolution of the fast solar wind streams in the two 
environments. 
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1. Introduction 

The solar wind is an inhomogeneous, anisotropic and compressible magnetized plasma 
where both velocity and magnetic fields fluctuate over a broad range of frequencies and scales, 
see e.g. the reviews of Tu & Marsch (1995) and Horbury & Tsurutani (2001). Fluctuations 
may originate either from the nonlinear interactions between large-scale streams (Coleman 
1966, 1968; Matthaeus et al. 1990) or by interacting Alfven waves produced close to the 
Sun and carried by the wind, (Belcher & Davis 1971; Dobrowolny et al. 1980; Leamon et 
al. 1998). Observations of the radial evolution of magnetic flclds in the inner heliosphere 
show the presence of fully developped turbulent spectra within a range of 10~^ — 10~^ Hz 
(Bavassano et al. 1982). 

The spectral index depends on the frequency range and on the distance from the 
Sun, varying from —1.2 to —1.7. Low-frequency measurements are performed at around 
10~^ — 10~^ Hz (Coleman 1968), while high frequency measurements sample the range closer 
to 10^2 - 10"^ Hz, (Bavassano et al. 1982; Leamon et al. 1998; Horbury & Balogh 2001). 
The spectral index tends to flatten closer to the Sun, indicating that turbulence is evolving 
in the solar wind. Phenomenological theory of hydrodynamic turbulence (Kolmogorov 1941) 
predicts a value of —5/3 for the spectral index, while the theory of Alfven wave-driven mag- 
netohydrodynamic (MHD) turbulence of Iroshnikov (1963) and Kraichnan (1965), predicts 
a slope of —3/2. Neither prediction takes into account the possbile influence of anisotropies 
and the presence of intermittency (Burlaga 1991, 1992; Marsch & Liu 1993; Carbone 1993; 
Feynman & Ruzmaikin 1994; Carbone, et al. 1995b; Horbury&; Balogh 1997; Ruzmaikin et 
al. 1995; Bruno et al. 2003; Hnat et al. 2003; Bershadskii & Sreenivasan 2004), in a systematic 
way. 

The presence of anisotropy makes it difficult to compare observed data with the two 
predictions, while the presence of intermittency tells us that the characteristics of the spec- 
trum are not sufficient to characterize the system: higher order statistics need to be taken 
into account. In particular, spectral indices alone are insufficient to discriminate amongst 
turbulence models. 

Anisotropy has been measured by various techniques involving the calculation of second 
order moments of the field either in the real or Fourier space, such as the variance matrix 
or the power spectra (Belcher & Davis 1971; Carbone et al. 1995). The eigenvector of the 
variance matrix corresponding to the minimum eigenvalue is usually known as the minimum 
variance direction. This direction is aligned with the large-scale mean field, indicating 
suppression of turbulence in that direction (Leamon et al. 1998; Bruno et al. 1999). Several 
MHD models incorporate at various levels the asymmetry of the spectral indices in the field- 
aligned (longitudinal) and transverse directions (Shebahn et al. 1983; Zank & Matthaeus 
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1992; Ng & Bhattacharjee 1996; Goldreich & Sridhar 1997; Matthaeus et al. 1998). Although 
this is a possible way to characterize anisotropy, second order longitudinal and transverse 
structure functions contain both anisotropic and isotropic contributions, as will be detailed 
later in Section 2. Those two contributions, always mixed, need proper treatment to be 
disentangled. A more systematic approach to analyse anisotropy is therefore important. 
Moreover, the relation between anisotropy and intermittency has not been investigated so 
far. Anisotropy and intermittency may also be important in the context of scattering of 
particles in the hehosphere (see e.g. Giacalone & Jokipii (1996)) 

Wc present in this paper a method for extracting in a systematic way, from the one- 
dimensional spacecraft data, information on the anisotropy and intermittency of the magnetic 
field fluctuations, and the interplay between them. We base our analysis on the behaviour 
of both diagonal and non-diagonal components of higher order structure functions. We have 
systematically compared isotropic and anisotropic fluctuations at different scales and for 
different magnetic correlation functions. We measure how fast isotropy is recovered at small 
scales, concerning both typical fluctuations of the order of the mean standard deviation, and 
highly intermittent events, affecting more the tails of the magnetic field probability density 
at all scales. We use Ulysses data of high speed streams at two different points along its orbit, 
at high and low latitudes, in order to assess the dependence on the large-scale properties of 
the small-scales anisotropic fluctuations, i.e. the issue of small-scales universality. 

The paper is organized as follows. In Section 2 we present the set of observables needed 
to have a systematic control on the isotropic and anisotropic ensembles. In Section 3 we 
present our data set and in Section 4 the main results for both the low and high latitudes 
data. Section 5, summarizes our flndings suggesting further possible investigations. 

2. Anisotropy and Structure function analysis 

Structure function decomposition into isotropic and anisotropic components has already 
been exploited with success in hydrodynamics, both for experimental and numerical data 
analysis (Arad et al. 1998, 1999; Kurien & Sreenivasan 2000; Biferale & Toschi 2001; Biferale 
& Vergassola 2001; Shen & Warhaft 2002b), see also Biferale & Procaccia (2004) for a recent 
review. In the latter case, the anisotropic contents of fully developed flow has been systemat- 
ically analyzed. Anisotropic fluctuations of the velocity fleld have been shown to be charac- 
terized by anomalous scaling, thus explaining the higher than predicted anisotropy found in 
the gradient statistics, an effect known as "smal-scales persistence of anisotropics" (Shen & 
Warhaft 2000, 2002; Biferale & Vergassola 2001). Some analytical results on the persistency 
of anisotropics for magnetic fields have also been obtained in the simplified case of passive 
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magnetic advection by stochastic velocity fields (Falkovich et al. 2001; Lanotte & Mazzino 
1999; Arad et al. 2000). The ideal way to assess the relative isotropic/anisotropic contents 
at all scales is to perform a decomposition of the correlation functions, of order 2 and higher, 
over a suitable eigenbasis with definite properties under the group of three-dimensional spa- 
tial rotations (the SO (3) group), corresponding to spherical harmonics decomposition for the 
simpler case of scalars. In principle, one needs to distinguish among different anisotropic 
contributions corresponding to the different projections on the whole set of eigenf unctions 
(Arad et al. 1999). Spacecraft data are inherently one-dimensional, therefore not directly 
suitable to be fed into an SO (3) analysis, which requires the whole field in a 3D volume to 
be systematically worked out. However, we shall show how it is possible to construct corre- 
lation functions of different orders that have null projection over the isotropic ensemble, i.e. 
with a leading contribution coming only by its leading anisotropic content, if any (Kurien & 
Sreenivasan 2000; Staicu et al. 2003; Jacob et al. 2004). 

Data analysis is based on a set of multi-scale correlation functions built upon differ- 
ent combinations of magnetic field components. The most general nth order correlation, 
Sai,...,a„{r), depending on single separation (r), is built from the n spatial increments of 
magnetic field components: 

^S...,.„(r-) = {5rB^,5rB^, ■ ■ ■ 5rB^„) (1) 

where 

6rB^ = B^{x + r)-B^{x) (2) 

is the difference between the values of component Ba at two different points a distance r 
away. Brackets (■) in (1) indicate the average over the locations x. Notice that in (1) we 
have assumed homogeneity but not isotropy, i.e. the correlation functions keep their explicit 
dependence on the full vector r. The correlation function (1) includes both isotropic and 
anisotropic contributions: 

SS' a (r) = S^^^'^Z {r) + S^^^'Z'^r). (3) 

In principle, different anisotropic contributions exist and one would need to further distin- 
guish among them. In this study we limit ourselves to disentangling the isotropic contribu- 
tions from the anisotropic, without entering the more subtle problem of separating out all 
the different anisotropies (the interest reader may consult Biferale & Procaccia (2004) for a 
detailed illustration on how to proceed in this direction). 

For n — 2 and ai — a2 in (1), we get the well known positively defined second order 
structure function, connected to the the energy spectrum Ea,a{k) — (|-BQ,(fc)p) via a Fourier 
transform. Another widely used form of (1) is the longitudinal structure function, obtained 
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by projecting all field increments along the separation versor, r: 5'£(r) = {{SrB • f)"). The 
general form of the tensor (1) for n = 2 in the case of a fully isotropic and parity invariant 
statistics, is given by the combination of the separation vector r and the only isotropic 
second order tensor, the unity matrix Sa,)3'- 

SSi7(r) = {SrBa,SrB^,r" = a(r)(5„,,„, + 6(r)r„,r„, (4) 

where a(r) and 6(r) are two scalar functions depending only on the amplitude r = \r\. 
Similarly, the expression for the fourth order isotropic tensors, Sai,'^-'^"a4{r), comprises three 
scalar functions, c{r),d{r), f{r): 

c(r)(5ai,a2(^a3,a4 + perm.) + d{r){5o,^^a2'razrai +perm.) (5) 

Analogous expressions hold for higher oder isotropic correlation functions. The key observa- 
tion is that by a suitable choice of the combination of indices ai, ■ ■ ■ , a„ and of the orientation 
r one may have the isotropic components vanish at any order, n in (3). From now on, let us 
fix the separation distance in the direction x so that r = (ra,, 0, 0). For the case n — 2, when 
7^ CK2, the resulting isotropic components vanish. We therefore have three different second 
order correlation functions that are purely anisotropic. When the order n of the correlation 
function is even it is enough to take an odd number of field increments in two different direc- 
tions to have a purely anisotropic observable. Therefore a possible set of purely anisotropic 
correlations have the form: 

SZi^^) = (^rM^r.Bj) {p + q = n) (6) 

with both p and q odd and such that p + q = n. The above nth order correlation has a 
vanishing isotropic component when the combinations of indices a — x and (5 — y,z are 
taken. 

Before presenting the results of our data analysis, let us briefly comment on the trans- 
lation from time series to spatial signals in our dataset. Of course, as it is the case for all 
spacecraft data, we only have access to the time evolution of the magnetic field along the 
orbit. We therefore cannot make an explicit evaluation of simultaneous field increments over 
space. Nevertheless, the advecting velocity speed is so high (see next section for a summary 
of the main physical relevant quantities) that in the range of frequencies we are interested in, 
it is possible to safely adopt the Taylor hypothesis and translate time increments into spatial 
increments. The "Taylor hypothesis" consists in supposing the 3D field as frozenly advected 
by the underlying large scale velocity field, Vq Frisch (1995). Field increments in the same 
spatial point at two times, t, f, are considered equal to the instantaneous field increments 
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over two spatial locations x and x-\-r with r — Vq {f — t). Therefore, for us, the direction 
r is fixed and given by the direction of the wind at the location of the spacecraft that is, 
within a few percent, the spacecraft-Sun direction. This direction, as said before, will be 
taken as our reference x axis. Spatial homogenity is translated via the Taylor hypothesis 
into temporal stationarity. 



Ulysses orbit samples the interplanetary plasma at distances varying approximately 
from 1 to 6 A.U, on a polar orbit. It is therefore possible to follow the evolution of plasma 
characteristics with distance and latitude. We use two different set of data: the first one 
was taken by Ulysses during 1992-1993, when the spacecraft was at about 20° heliographic 
latitude and 5 AU distance from the Sun. The second was taken at the end of 1994, with 
Ulysses above the South Pole, at about 80° latitude and a distance of about 2 AU from 
the Sun. This latest dataset has just recently been made available to the community by 
the Ulysses team. Solar activity was, during the 92-93 period, dechning, after the 1990 
maximum. In 1994, the cycle was approaching the minimum of 1996. Each daily dataset 
provides the magnitude of all three components of the interplanetary magnetic field, taken 
at the rate of 1 or two seconds by the Vector Helium Magnetometer on board (Balogh et 
al., 1992). In Table 1 we report, for the two datascts, the interval of time considered, the 
heliographic latitudes spanned, distance from the Sun, average speed of the wind and average 
magnetic field intensity. 

We pre-process data in order to clean spikes due to instrumental problems or to large 
shocks. This is made by excluding those data where the jump in the magnetic field between 
two consecutive data points (usually 1 second apart) is larger than a threshold, AB, of the 
order of the mean large scale magnetic field. A fraction of datapoints as small as 10~^, 
is discarded this way. As a result, we can access magnetic field fiuctuation on a range of 



3. Ulysses Dataset 



Table 1. Low latitudes and Polar datasets. 



Dataset 



Days 



(HGL) 



Lat 



(A.U.) 



Dist 



Speed {B) 
(Km s-i) (nT) 



Low lat 92/209-93/137 -15 to -30 5.3 to 4.7 750 0.47 
Polar 94/245-265 -79.7 to -80.2 2.37 to 2.23 760 1.3 



- 7- 



frequencies of almost six decades. 

In Table 2 we detail the total number of datapoints in the dataset, the number of datapoints 
discarded, iVexci, the fraction of the latter to the total, the threshold on the maximum jump 
between magnetic field for consecutive datapoints, the average field intensity for the whole 
dataset. 



The alternating pattern of slow and fast wind is shown in Fig. 1, spanning a ten month 
period, from day 209, 1992, to day 137, 1993. Within this period, wc selected those sequences, 
of about five days each, when spacecraft is embedded in the trailing edges of high speed 
streams and velocity is above 650 Km/s. The days selected arc, in 1992, 209-214, 235-241, 
259-263, 337-342 and, in 1993, 28-34, 53-57, 81-85, 108-113, 133-137. They are highlited in 
Fig. 1 within vertical lines. 



Twenty-one consecutive days around the maximum latitude reached at pcrihclium, dur- 
ing the fast latitude scan of 1994, arc selected. Differently from the previous dataset, only 
the fast component of the wind is present. Table 1 lists latitude range, distance, average 
speed and average magnetic field for this dataset as well. 



3.1. 



Low latitude dataset 



3.2. High latitude dataset 



Table 2. Data selection 



Dataset 



N 



A^excl A^excl/A^ (B) 

{nT) {nT) 



Low latitude 3 915 792 



78 



2.0 E-5 



0.5 



0.47 



High latitude 1 476 051 36 2.4 E-5 1.2 1.30 
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Fig. 1. — Plasma velocity sampled by Ulysses spacecraft between day 209 (July 27) 1992 
and day 137 (May 17) 1993. Spacecraft was between —15° to —30° heliographic latitude, 
approaching the Sun at a distance varying from 5.3 to 4.7 AU (see Table 1). Vertical lines 
highlight selected intervals in the trailing edges of high-speed streams. 
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4. Results 

4.1. Equatorial data 

We want to first test the consistency between the disjoint sets making up the low latitude 
dataset of Fig. 1. The second order longitudinal structure functions, calculated for each of 
those intervals of contiguous data, are shown in Fig. 2. They are consistent with each other 
over more than 5 decades, from 1 to 10^ Hz in the spacecraft frame, which translates, with 
a mean plasma velocity of 750 Km/s, into a range of 7.5 • 10~^ Mm to 7.5 • 10^ Mm. Some 
intervals have a more intense signal than others do. 

The anisotropic component Sxz shown in the inset of the same figure, displays a similar 
behavior. Wc conclude that data from different intervals are commensurable and combine 
them together to obtain more stable statistical results. We shall refer to the combined set 
as the "low-latitude" dataset without further distinction. 

Let us now compare the undecomposed second order structure functions with its anisotropic 
content. In Fig. 3 we plot the longitudinal structure functions of second order, S^}x{tx) and 
the two transverse structure functions in the directions perpendicular to the x axis, Syy {rx) 
and Szz [tx)- All these functions have both isotropic and anisotropic contriution: 

= S'i^ri.Tx) + 5S'^"^^°(r.) (7) 

The two purely anisotropic second order structure functions S^3y{rx) and S^){rx), are plotted 
in the same figure. A few comments are in order. First, we notice that the anisotropic corre- 
lations have a smaller amplitude with respect to the full correlation functions. This suggests 
that the isotropic contribution in the decomposition (3) is dominant. Moreover, we see that 
the anisotropic curves decay slightly faster than the full correlation by decreasing the scale. 
In other words, isotropic fluctuations become more leading going to small scales, but very 
slowly. This is consistent with the recovery-of-isotropy assumption often advocated in many 
phenomenological theory of hydrodynamic turbulence and magnetized plasma. However, in 
order to assess more precisely this issue, it is important to control higher order statistical 
objects, i.e. the whole shape of the probability density distribution, at all scales. In the inset 
of the same Fig. 3 we show the same comparison between longitudinal, Si^!xx{rx), transverse, 
Saaaa{rx) (with a = y,z) and purely anisotropic correlations oi fourth order (see caption in 
the figure). Now the situation is quite different. First, the intensity of some purely anisotropic 
components are much closer to those with mixed isotropic and anisotropic contributions, i.e. 
the longitudinal and transverse structure functions. Second, the decay rate as a function of 
the scale is almost the same: no recovery of isotropy is any more detected for fluctuations of 
this order. A similar, even more pronounced, trend is observed for sixth order quantities (not 
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[Mm] 

Fig. 2. — Second order longitudinal structure function, S^J(rx), for each interval comprising 
the low latitude dataset (see Fig. 1), as a function of the separation r^- In the inset, the 
second order purely anisotropic structure function, Sxz {tx) ■ 
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shown). The persistence of strong anisotropics at high frequencies (small scales) cast some 
caveat on measurements of quantities which do not properly disentangle the isotropic from 
the anisotropic components. As it will be shown later for the case of high-latitudes data, 
anisotropic components have strong variations in intensity depending on the position on the 
solar orbit. Both latitude and distance from the Sun influence the amount of anisotropy. 
As a result, undecomposed quantities which are influenced by both isotropic and anisotropic 
fluctuations are expected to be non-universal, the anisotropic content being dependent on 
the spacecraft position and latitude. This must hold for the spectrum and even more for 
higher order structure functions. 



Anisotropic fluctuations are not the unique source of complexity in solar wind data. It 
is well known that both magnetic and velocity flelds are strongly intermittent, i.e. their 
statistical properties at different scales cannot be simply superimposed by rescaling. This 
implies the existence of anomalous scaling laws in the structure functions and "fat tails" 
in the PDFs of held increments Prisch (1995). Here we want to address this issue for the 
anisotropic sectors. The main conclusion will be that anisotropic correlations also show 
anomalous scaling, their PDFs becoming more and more non-Gaussian at small scales. In 
Fig. 4 we show the Kurtosis of both the longitudinal and transverse structure functions, 
i.e. the ratio between fourth order moments and square of the second order moments of 
longitudinal and transverse increments: 



A Gaussian variable would have a Kurtosis of 3, independent on the scale while all three 
curves grow at small scales. We stress once more here that these quantities probe both 
the iso and anisotropic physics. Therefore the scaling properties are certainly affected by 
the superposition of different contributions. In the previous section we have shown that 
the isotropic sector is never sub-leading. We may therefore consider the above result as a 
conflrmation that the isotropic fluctuations are indeed strongly intermittent. 

Similarly, to investigate intcrmittency in the anisotropic sector, it is useful to deflne a 
purely anisotropic Kurtosis, by taking the adimensional ratios of fourth order and second 
order anisotropic correlation functions: 



4.2. Intermittency 





(8) 



{Saa (^x))^ 
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[Mm] 

Fig. 3. — Second order longitudinal, transverse and purely anisotropic structure functions. 
Low latitude dataset. The upper three curves show the longitudinal and tranverse structure 
functions: solid line — S^^x] empty cirles o S'i^y] filled circles • S^zz ■ Errorbars are superim- 

(2) 

posed on — Sxx- Errors are evaluated as the standard deviation of the individual intervals 
comprising the whole dataset. The lower curves show the purely anisotropic structure func- 

f2) f2) f2) 

tions: Sxy, ▲ filled triangles; SxJ ■, A empty triangles; Syz , □ empty squares. Errorbars are 
superimposed on A S'i) . Inset: fourth order structure functions, longitudinal, transverse and 
purely anisotropic. Solid line, — S'^xxx', empty circles o S^yy] filled circles • S^zzz- Purely 
anisotropic structure functions are: Sx%y, A filled triangles; Si%z, A empty triangles; Sy^zz, 
□ empty squares. 
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where a, P are chosen so that contributions from the isotropic sector in both the numerator 
and the denominator vanish. The anisotropic components of the kurtosis (9) are shown in 
the inset of Fig. 4. Functions are increasing towards small scales, with slopes of t^^™*" = 
-0.6 ± 0.2, xT''" = -0-8 ± 0.2, xT''° = -0-45 ± 0.2 for the xy, xz and yz components, 
respectively (see Table 3). This is the first clear indication, to our knowledge, that anisotropic 
fluctuations in the solar plasma are strongly intermittent. Similar trends are observed for 
generalized Kurtosis of sixth order (not shown): 



(aniso) 



There, our best estimate for the exponents is Xg = —1.2 ± 0.3, xy component, and 
Xe""**"'* = — 1.5±0.3, 0:2; component. 

Let us here remark that the quantity in (9) is not constructed from ratios of 4th and 
2nd order moments of the same observable, i.e. it is not, rigorously speaking, the kurtosis of 
a stochastic variable. Nevertheless, it is a good probe of the relative intensity of 4th versus 
2nd order anisotropic moments, the best that can be done with a one- dimensional set of 
data. 

A power law fit of the numerator and denominator of (9,10) can be used to directly 
measure the scaling exponents of the second order, 

S^S^r^)-ri^ , (11) 
and higher order anisotropic correlation functions, 

/A\ , , J aniso) ( a\ , , J aniso) 

tUir^) ^ ri^ ^i%,(r.) ~ r^^ (12) 

with, as customary now, a, (3 arc chosen in such a way that only purely anisotropic quantities 
are returned. We found (^^"™'^°) = o.75 ± 0.1 for the xy component, (^^"""°) = o.95 ± 0.1 for 
the xz component, and (^^"™^°) = 0.75 ± 0.1 for the yz component, see Table 3. Values for 
the fourth and sixth orders (^|"™*°) ^^""-^s") j^g^y g^jg^ read out from the same table. 
Errorbars are estimated from the change of slope in the range of scales from 10 to 10^ Mm. 
Missing entries in the table indicate that the scaling properties were not well defined within 
that range. 

The above results show that anisotropic fiuctuations, although they never become the 
leading ones are still important at small scales. Order by order, the undecomposed correlation 
function is more intense than any anisotropic projection. This can be visualized, for the 4th 
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Fig. 4. — Kurtosis (8) of longituinal and transverse magnetic field fluctuations. Solid line — 
K^x\rri.)] empty circles o Klpirr,.)] filled circles • K^^\rx)- Straight lines represent linear fit 
to the central portion of the curves with slopes of —0.31 for the longitudinal component, and 
—0.38 for the two transverse ones. The horizontal line corresponds to the Gaussian value 
of 3, attained only at large scales. In the inset: purely anisotropic kurtosis (9). K'^y{rx)-, A 
filled triangles; K^){rx), A empty triangles; K^z{tx)i D empty squares. Straight lines have 
slopes of -)(^'^^° = — 0.8 and -)(^'^^° — —0.6 in the xz and xy components, respectively. Low 
latitude dataset. 
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and 6th orders, by plotting the ratio between the undecomposed object and one anisotropic 
projection: 

'^xz {^x) - (4) , '^xz [r-x) - -(Q) — U^J 

iJxxxxyx) iJxxxxxxyx) 

These quantities never increase at small scales, indicating that isotropic contribution in the 
denominator is leading with respect to the anisotropic, see Fig. 5. Another quantity that 
can be used to characterize the relative weight of anisotropic to isotropic fluctuations, may 
be built from a nth order anisotropic moment and the n/2 power of a 2nd order isotropic 
moment (Shen & Warhaft 2002b; Biferale & Vergassola 2001). For example, in our geometry, 
one possible choice would be: 

0(4) / \ q(6) / N 

F^f (rx) = Zr''^ F(f (r.) = B^^IfM (14) 

where the numerator is a purely anisotropic nth order quantity while the denominator is the 
2d order longitudinal structure function, raised to the n/2 power. Clearly all quantities in 
(13) and (14) would be vanishing in a perfect isotropic ensemble. The difference between the 
two deflnitions (13) and (14) for F and G, lies in the normalizing function in the denomina- 
tor. In the flrst case, G, the normalization is through a correlation of the same order of the 
numerator while in the second case, F, is via a second order correlation raised to the appro- 
priate power. Their amplitude as a function of rx can be taken as a measure of the change in 
the anisotropic content as a function of scale. The defintion (14), on the other hand, mixes 
correlation of different orders, thus including their possible different intermittent corrections 
(Biferale & Vergassola 2001). In Fig. 5 we also show the behavior of Fx^\rx) for n = 4,6. 
Again, there is a clear indication of the presence of important anisotropic contributions, 
particularly at small scales. 



4.3. Probability density functions 

Before concluding this section we want to re-discuss some of the previous results from the 
point of view of the probability density functions (PDFs). Anisotropics may be highlighted 
at the level of the PDF by looking at the antisymmetric part of the distribution of field 
increments at different scales. Let us define the PDF, P(Xq,^), of the dimensionless magnetic 
field increments at scale r^,: 

X^fsirx) = ^^^-^"'^^^-^/^ (15) 

{^r^Bx5r^Bx) 

In order to make the stochastic variable dimensionless we have normalized it with the longi- 
tudinal second order structure functions at that scale. With a suitable choice of the indices 
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Fig. 5. — Generalized flatness (j^I^^ ^^"^ ^a^^^^ of order 4 and 6, Eqs. (13) and (14) for 
components xz and xy. Low latitude dataset. 



q;/3, all odd moments of Xapivx) would be zero in a perfectly isotropic ensemble. This is the 
case when a — x and (3 — y,z. We may now define the antisymmetric part of P(Xq^) as 



(16) 



and notice that it would vanish in a symmetric isotropic ensemble. 

Ar{Xai3) gives us a direct measurement of the anisotropy as the imbalance in the prob- 
ability of having oppositely directed fluctuations at that scale. In Fig. 6 we show the anti- 
symmetric part of the PDF, Ar{Xai3{rx)) for a = x and (3 = z for three different separations 
fx- The increasingly fat tails as one goes to smaller scales, reflects the non-gaussianity of 
P{Xxz{tx))-, which becomes more enhanced at small scales. In order to assess the relative 
weight of the antisymmetric versus the symmetric fluctuations, we deflne the normalised 
antisymmetric part of P(Xa^): 



This quantity also vanish in a symmetric isotropic ensemble, approaching the value one in 
the limit case of strong anisotropy, P(Xq,^) ^ P{—Xaf3). In the inset of the same flgure, 
R[Xxz) is shown. The fact that at large separations R{Xxz) is close to one, means that large 
events are progressively more anisotropic as they grow in intensity, a possible signature of 
the large scale structures in the plasma. For small separations, the system is indeed globally 
more isotropic, although small scale anisotropy never vanish and survives at a significant 
level of 10% for all intensities. 



We discuss here anisotropy and intermittency detected in the polar region by Ulysses. 
This allows us to address the "universality" of anisotropy, i.e. quantifying to which extent 
intensities of anisotrpic fiuctuations and their scaling properties are dependent/independent 
on the mean large scale structure on the magnetized plasma. There are two effects which 
might influence the relative anisotropy of the turbulence in the polar and equatorial regions. 
In the polar regions, the amplitude of turbulence relative to the mean fleld is stronger, while 
the effects of solar rotation, which tend to bend the interplanetary magnetic fleld into a 
spiral, are negligible. In the equatorial high speed streams, the average magnetic fleld is 
bent into the Parker (spiral) direction, so that there are two main axes which may influence 
the evolution of the fluctuations, the radial and the mean fleld directions. The mean fleld 
direction coincides with the radial direction for polar flows while it is perpendicular to it, 
close to the y direction, for the low-latitude data arlund 5 A.U. 




P{X^p{rx)) - P{-X^p{Tx)) 



(17) 




4.4. High-latitudes data 
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Fig. 6. — Antisymmetric part of the PDF of X,j:z{rx), Ar{—Xxz) Eq. (16), for three different 
spatial separations r. Sohd hne: r = 12 Mm, dotted hne: r = 192 Mm, dot-dashed Une: 
r = 3072 Mm. Inset: the normahsed antisymmetric part of the PDF, Rx{—Xxz), Eq. (17), 
for the same set of rx- 
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Let us first present results on the overall relative importance of anisotropic fluctuation 
with respect to the undecomposed ones. In Fig. (7) we show the same as in Fig. (3) but 
for polar data. Purely anisotropic structure functions have a much lower intensity (one 
order of magnitude less) with respect to the longitudinal and transverse structure functions 
both for the second order (body of the figure) and for the fourth order (inset). Indeed, for 
higher order moments, 6 and higher, the statistical fluctuations combined with the very low 
intensity of the anisotropic signal do not allow to have stable results even with the whole 
statistic of 21 consecutive days we analyzed. We conclude therefore that the anisotropy 
content at this latitude is much lower than in the low latitude dataset. One could argue that 
at this latitudes averaging over long periods may hide important physical phenomena which 
appear on a shorter time window. Therefore, we also selected periods of 2-3 consecutive days 
when the anisotropic signal looked more stable and intense. The anisotropic content in those 
events is sligthly more important and allow to make a quantitative estimate of its scaling 
properties, but do not differ qualitatively. In Table 3, data with an asterisk * indicate that 
scaling exponents are evaluated on the smaller dataset. 

In Fig. (8) we show the same as Fig. (4), for the polar data set. We show the Kurtosis of 
longitudinal and transverse magnetic field fluctuations toghether with the Kurtosis for purely 
anisotropic correlation functions (9). Comparing the scaling behaviours of all the statistical 
indicators considered, summarized in Table 3, we have a qualitative agreement between the 
polar and the "Equatorial" data set. If confirmed by other measurements, and/or with 
higher statistical data sets, this would be a nice indication of "universality" in the small 
scales fluctuations of the solar wind plasma. Overall intensities of isotropic and anisotropic 
contents are of course dependent on the distance and latitude, while their variation with 
scale/frequency look more stable. 



5. Conclusions 

Our main finding is the detection of strong anisotropic fluctuations in the equatorial part 
of the orbit. Here, the anisotropic contents of fourth order correlation function is roughly of 
the same order of its isotropic part, at all scales, indicating that small scale isotropy is not 
recovered. Moreover, a high degree of intermittency is measured in the purely anisotropic 
fluctuations. In the polar region, anisotropics are smaller and highly fluctuating in time, but 
with a spatial dependencies compatible, within statistical errors, with the one observed at 
low latitudes. This would indicate some universal features of anisotropic solar fluctuations 
independently of the latitude, at least for what concerns their scaling properties. Our results 
point toward a crucial role played by anisotropic fluctuations in the small scales statistics. 
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[Mm] 

Fig. 7. — Same as in Fig. (3). Polar dataset. Second order longitudinal, transverse and purely 
anisotropic structure functions. The upper three curves show the longitudinal and tranverse 
structure functions: solid line — S^x , empty cirles o S^y and filled circles • S^^z ■ The 

(21 

lower curve with triangles show the largest purely anisotropic structure functions: Sxy ■, A 

(2) 

filled triangles; Syz , □ empty squares. Inset: fourth order structure functions, longitudinal, 
transverse and purely anisotropic. Solid line, — S'^xxx, empty circles o Sy%y, filled circles • 
Sifzz- Purely anisotropic structure functions are: Si%y, A filled triangles; Si%z, A empty 
triangles; Sy'^zz, D empty squares. 



- 21 - 
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[Mm] 

Fig. 8. — Same as in Fig. (4). Polar dataset. Kurtosis (8) of the longituinal and transverse 
magnetic field fluctuations. Solid line — K^\rx)-, empty circles o K'^\rx) and filled circles • 
K^f' (fx) ■ The dot-dashed hne shows the constant level of three for the kurtosis of a Gaussian 
variable. In the inset: purely anisotropic kurtosis (9) of component A Kx'^(rx)- 
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Models where higher order statistic is also taken into account, providing estimates for the 
scaling exponents of higher order anisotropic structure functions, will be important to a 
deeper understanding of solar wind turbulence. 

Before concluding, let us go back to the issue of distinguishing different anisotropic 
fluctuations. As one learns from the theory of group of rotation in three dimensions, there 
is not a unique "anisotropic" sector, rather different anisotropic properties are described by 
projection on the eigenfunctions with different total angular momentum, j, and projections 
of the total angular momentum on a given axis, m, Arad et al. (1999). As mentioned in 
the introduction, the exact decomposition in different anisotropic sectors is possbile only if 
using numerical data, providing access to the whole magnetic field in the 3D space. Here 
we have described the procedure that should be adopted for one-dimensional strings of data, 
extracting the "whole" anisotropic components out of the experimental data. This implies 
that all the estimate of scahng properties as here reported may well be affected by out-of- 
control contributions from different anisotropic sectors. The hope is that out of all anisotropic 
sectors, only the leading one is dominating the statistics at small scales. This hypothesis, 
which implies a hierachy between the scaling exponents in different sector has been verified 
on direct numerical simulations of turbulent flows Biferale & Toschi (2001); Biferale et al. 
(2002) and on analytical calculation for passive magnetic fields Lanotte & Mazzino (1999); 
Arad et al. (2000), but remains an open question for active magnetic fields. 

We are thankful to B. Bavassano, R. Bruno, A. Lanotte and F. Toschi for friutful 
discussions. We acknowledge support from EU under the grant "Nonideal Turbulence" 
HPRN-CT-2000-0162. 
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Table 3. Universality. Scaling exponents, Eqs. (9,11,12) 



/^aniso 
S2 


/^aniso 
S4 


/-aniso 
S6 


^. aniso 
X4 


^.aniso 
A6 


Low-lat 0.75±0.15(a;y) 
0.95 ±0.10 {xz) 
0.75 ±0.10 {yz) 


0.8 ± 0.3 {xxxy^ ) 
1.0 ±0.15 {xzzz) 
1.0 ± 0.25 [yzzz) 


1.2 ± 0.4 {xxxyyy) 
1.2 ± Q.2{xxxzzz) 

2 ± 1 (yyyzzz) 


-0.6 ± 0.2 (xy) 
-0.8 ±0.2 (xz) 
-0.45 ± 0.2 (yz) 


-1.2 ± 0.3 (xy 
-1.5 ±0.3 (xz 


Hi-lat 0.75 ± 0.15 (xy) 
0.75 ±0.15 {yz) 


0.8 ± 0.2 {xxxy"^) 
0.8 ± 0.3 {xzzz) 
1.0 ± 0.3 {yzzz) 


1.1 ± 0.3 (xxxyyy) 
1.1 ± 0.3 (xxxzzz* ) 
1.5 ± 0.4 (yyyzzz* ) 


-0.6 ±0.2 (yxt) 


-1.1 ±0.3 (xy 



*Exponents evaluated on the short polar dataset, see Sect. 4.4. 

^First component provides a third order correlation rather than the second one, as in all other cas 
considered. 



